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Abstract 

We study the origin of baryonic matter via leptogenesis in realistic SO(10) models, in particular, 
in a new lopsided mass matrix model introduced recently by three of the authors. By introducing 
simple CP-violating phases in the mass matrix of the right-handed neutrinos, the model generates 
sufficient baryon asymmetry without fine-tuning. We compare this result with other realistic 
SO(10) models. 
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I. INTRODUCTION 



One of the most fundamental questions in modern cosmology is where the baryon number 
asymmetry in today's Universe comes from. It is an attractive idea to assume that it 
arises dynamically in a symmetric big bang model through C and CP violating processes 
out of thermal equilibrium. Of the several proposals that use this scenario, leptogenesis 

has emerged as one of the most interesting scenarios for baryon generation. Here the 
baryon number is produced through the so-called sphaleron processes from a residual lepton 
number density left over from an early stage of the universe through lepton-number violating 
decays of right handed neutrinos. Such connection is extremely interesting since the same 
heavy right-handed neutrinos are also important ingredients for understanding the small 
left-handed neutrino masses through the so-called see-saw mechanism |2|. 

Recently much study has been made in the literature about the feasibility of such a 
scenario |3|. Although one can discuss leptogenesis and neutrino mass texture in models 
without introducing quark degrees of freedom, the most interesting theoretical frameworks 
involve quark-lepton unification, possibly under supersymmetric formalism. In addition to 
connecting the origin of matter to the big picture of particle physics, the grand unified 
theory (GUT) models build in more constraints and have better predictive power For 
various reasons, SUSY SO(10) has been a favored framework to unify physics beyond the 
standard model (for a recent review of SO(10) models, see 0). Depending upon which 
set of Higgs multiplets is chosen to break the GUT group and electroweak symmetry, two 
classes of SO(10) models are most studied in the literature: One uses IOh, 126^, 126# and 
120 H 0,0,113, and the other uses 10 H , 16 H , 16h and 45# 0,QE3,E3- While most of 
these models are quite successful in fitting and predicting the known experimental masses 
and mixing angles of leptons and quarks, they predict very different values for the poorly- 
known neutrino mixing angle 8 13 . Quite interestingly, they also paint different pictures for 
leptogenesis. 

In this paper, we are mostly interested in the leptogenesis in a lopsided SO(10) mass 
matrix model proposed recently by three of us (X.J, Y.L. and R. N. M.) |l3|. The model is 
a modification of the lopsided model originally proposed by Albright, Babu and Barr 
In the modified version, the right-handed neutrino mass matrix has a simple diagonal form. 
The large solar mixing angle is mainly generated from the neutrino Dirac mass matrices. 
The prediction of this model for sin #13 is much larger than the original model. Here we will 
show that the leptogenesis picture is also strikingly different. In particular, enough baryon 
symmetry can be produced without requiring the two heavy right-handed neutrino masses, 
Mi and M 2 , to be quasi-degenerate. 

The presentation of the paper is as follows. In Section II, we review the new lopsided 
SO (10) model. In Section III, we consider in detail leptogenesis in this model. In Section IV, 
we make a comparison of the leptogenesis with other SO(10) models, emphasizing similarities 
and differences. We present our conclusions in Section V. 



II. A NEW LOPSIDED SO(10) MODEL 

A new SUSY SO (10) GUT model with lopsided mass matrices was introduced in [1J: 
In this section, we briefly describe its physics. The model was obtained by modifying the 
right-handed neutrino, up-quark, and left-handed neutrino mass matrices of the original 
lopsided model of Albright and Barr 2 3; which we will describe briefly in Section IV. 
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The lopsidedness built within the Yukawa couplings between the second and third families 
generates, among other interesting physical consequences, the large atmospheric-neutrino 
mixing angle #23 while keeping V c b in the Cabbibo-Kobayashi-Moskawa (CKM) matrix small. 
The modified lopsided model assumes that the right-handed neutrino Majorana mass matrix 
M R has a simple diagonal structure, and introduces additional off-diagonal couplings in the 
upper-type-quark and neutrino Dirac mass matrices to generate 1-2 (solar angle) rotation. 
All the fermion masses and mixing angles can be fitted well in the new model. The mixing 
angle #13, however, is close to the upper limit from the CHOOZ experiment and therefore 
definitely within the range of next generation reactor experiments. 

We use the convention that Yukawa couplings in the lagrangian appear as 

£ = —QiHYijdjij + ... . (1) 

Then the fermion matrix matrices are = vY^ with v = 174 GeV. Through couplings 
with a set of Higgs multiplets 10#, 16#, 16# and 45#, the up-type-quark, down-type- quark, 
charged-lepton and neutrino Dirac and Majorana mass matrices in the model of Ref. 
have the following forms, 



M u = uj \Au , M VD = I u I At 
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M VR = b \A R . (2) 
V° V 

From the above, we get the Majorana mass matrix of left-handed neutrinos from the see-saw 
formula, m VL = -M Ud M~^M^ d , 

(rf I 'a + (k + p) 2 (k + p)uj t)(k - p) /a + (k + p)\ 
(k + P )uj uj 2 uj M 2 /A R . (3) 

ij(k- p)/a + (K + p) uj l + («- p) 2 /a + uj 2 /bj 

The parameter a is of order one, signaling the lopsidedness between the second and third 
families in Md and Mi. This feature leads to a large left-handed neutrino mixing in the PMNS 
matrix and a small left-handed quark mixing shown in the CKM matrix. The parameter e 
is one order-of-magnitude smaller than a and generates the hierarchy between the second 
and third families. In extending to the first family, S and 5' were introduced into the M& 
and Mi. Usually a large rotation in the 1-2 sector is in conflict with the hierarchial property 
of the quark masses. However, in the above texture, the 1-2 rotation angle from M u will 
be combined with the 1-2 rotation from to obtain the Cabibbo angle 8 C , and thus the 
constraint from the up-type quark spectrum is avoided. The first two families in the M u 
and M Vh are not coupled to each other directly but through couplings with the third family. 
The rotations in 1-2 sector generated for left-handed up-type quarks and neutrinos are 
proportional to the ratios 7 = (k — p/3)/u and 7' = (k + p)/uj, respectively. 

The procedure to fit various parameters to experimental data is as follows. First, the up- 
type quark and lepton spectra and the parameters in the CKM matrix are used to determine 
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10 parameters: a, e, 5, 5', 0, u>, 7, 77, Au and Ad- The best fit yields a and e approximately 
the same as those in the original lopsided model, and thus the successful prediction for the 
mass ratios m°/m° and m®/m b are kept. The down-type quark mass spectrum comes out as 
predictions. The present model constrains the neutrino mass spectrum as hierarchial. The 
mass difference Am 2 12 is used to fix the right-handed neutrino mass scale A R . 

We summarize our input and detailed fits as follows. For CKM matrix elements, we take 
\V US \ = 0.224, \V ub \ = 0.0037, \V cb \ = 0.042, and 5 CP = 60° as inputs at the electroweak 
scale. With a running factor of 0.8853 for \V u b\, and \V c b\ taken into account, we have 
\V® b \ = 0.0033 and \V® b \ = 0.037 at the GUT scale. For charged lepton masses and up quark 
masses, we take the values at the GUT scale corresponding to tan/5 = 10 from Ref. [l6^ . 
For neutrino oscillation data, we take the solar-neutrino angle to be # S oiar = 32.5° and mass 
square differences as Am^ 12 = 7.9 x 10~ 5 eV 2 and AmJ 23 = 2.4 x 10~ 3 eV 2 . The result for 
the 12 fitted parameters is 

a = 1.83 , e = 0.1446 , 5 = 0.01 , 

8' = 0.014 , = 27.9° , 77 = 1.02 x 10~ 5 , 

u = -0.0466 , p = 0.0092 , k = 0.0191 , 

M v = 82.2 GeV , M D = 583.5 MeV , A R = 1.85 x 10 13 GeV . (4) 

There is a combined constraint on a and b, and thus the right-handed Majorana mass 
spectrum is not well determined. As examples, if a = b, a = —2.039 x 10 -3 ; and if a — 1, 
b = -1.951 x 10~ 3 . 

In our previous paper, we have taken a = b and real. The results for the down-type quark 
masses and right-handed Majorana neutrino masses are as follows, 

m° d = 1.08 MeV , m° s = 25.97 MeV , m° b = 1.242 GeV , 

Mi = 3.77 x 10 10 GeV , M 2 = 3.77 x 10 10 GeV , M 3 = 1.85 x 10 13 GeV . (5) 

The predictions for the mixing angles in the PMNS matrix are, 

sin 2 fl atm = 0.49 , sin 2 29 13 = 0.074 . (6) 

If one releases the best-fit value of Am 2 12 and Am 2 23 and impose only the 3cr constraint 
as 7.1 x 10~ 5 eV 2 < Am 2 12 < 8.9 x 10~ 5 eV 2 and 1.4 x 10" 3 eV 2 < Am 2 23 < 3.3 x 10" 3 eV 2 , 
one would obtain, 0.44 < sin 2 # atm < 0.52 which is well within the la limit, and 0.055 < 
sin 2 2#i3 — 0.110 which, as a whole region, lies in the scope of the next generation of reactor 
experiments. 



III. LEPTOGENESIS 



The baryon asymmetry of the universe is customarily defined as the ratio of the baryon 
density and the photon density after recombination, and has been measured to very good 
precision from the WMAP experiment [l7| : 

t]b = — — = (6.1 i 0.2) x 10 ^. (7) 
Interestingly, the big-bang nucleosynthesis is completely consistent with this determination. 



4 



To produce this asymmetry through leptogenesis, several considerations have to be ad- 
dressed. First, what is the number of right-handed neutrinos decaying out of thermal equi- 
librium? The answer to this question is in principle depends on the thermal history of the 
right-handed neutrinos. In our model, it turns out that this dependence is rather weak 
because of the strong washout. Second, what is the lepton density generated from the right- 
handed neutrino decay? This, of course, is related to the CP asymmetry of the decay which 
depends on complex phases in the Yukawa interactions. Third, some of the generated lepton 
density gets washed out by inverse-decay processes and scattering. This effect can be rather 
important, particularly in the so-called strong washout region. Finally, one must calculate 
the percentage of lepton number density converted into the baryon number density through 
the electroweak sphaleron process. The answers to some of the questions are less model- 
dependent and are standard in the literature [18| . Here we focus on the parts depending on 
particular models for the right-handed neutrinos. 

The density of leptons from right-handed neutrino decays is 

3C(3)^T 3 A 

UL = 4^ Z^' ( 8 ) 

i=l 

where the first factor is the thermal density of a relativistic fermion with g^ = 2 and the 
sum is over the number of right-handed neutrinos. The q is the decay CP asymmetry of 
the i-th right-handed neutrino; Kj is the corresponding efficiency factor, taking into account 
the fraction of out-of-equilibrium decays and the washout effect. Both factors depend on 
the effective mass defined as 

Mi v ; 

where M' denotes the Dirac neutrino mass in the basis in which the right-handed neutrino 
matrix is diagonal, M' v = M U U*. If U diagonalizes M VR , then 

(M x \ 

M UR = UMU T = U M 2 \U T , (10) 

V o o m 3 / 

and the right-handed neutrino mass eigenstates are \% = ^fiXfi where Xf is the family 
basis. Reversing this relation, one has Xf = Yli U%Xi- Thus if the right-handed field Xf is 
replaced by Xi i n the Yukawa coupling, the Yukawa matrix is multiplied on the right by U* . 
Hence the above mass relation follows. 

The lepton number is converted into the baryon number through the B — L conserving 
electroweak sphaleron effect |l9j . 

8N G + AN H 

nB = -22N G + 13N H nL ' (U) 

where N G = 3 is the number of fermion families and iV# is the number of Higgs doublets. 
In the standard model N = 3 and Nh = 1. 

The photon density can be calculated from the entropy density s = -^g^T 3 , where g* 
is the effective number of degrees of freedom, through the relation 

4vr 2 / 21 \ 7T 2 , s 

S = 330 ( 2+ llj 2CM^ ^ 
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where the second factor takes into account the neutrino contribution. Ignoring the lightest 
right-handed neutrino contribution, g* is 106.75 in the standard model. 

The final ratio of baryon to photon number density through leptogenesis is 

' i i 

Now we turn to the decay asymmetry and efficiency factors. 



A. The CP Asymmetry from Right-Handed Neutrino Decay 

The right-handed neutrinos are assumed to be CP eigenstates in the absence of the 
Yukawa type of weak interactions. In the presence of the interactions, they can decay into 
both left-handed leptons (neutrino and charged leptons) plus Higgs bosons and right-handed 
antileptons plus Higgs bosons. In the leading order, the decay rate is 

r< = ^(Y*Y%Mi , (14) 

where again, Y' = YU* is the Yukawa matrix in the basis where the right-handed neutrinos 
are in mass eigenstates. 

At next-to-leading order, the decay rates into leptons and antileptons are different due 
to the complex phases in the Yukawa couplings. The decay CP asymmetry is defined as 

r(iv t -> ijH) - r(iv t -> jjgt) 

In one-loop approximation, one finds, 

" ^jr/KW (y*y% ' (16) 



(17) 



where the decay function is given by |2 

F(x) = y/x 



1 , s 1 +x 

+ 1- 1 + xln 



1 — X X 



In the limit of large x, this become —3/2y/x. The first term in F is singular when two right- 
handed neutrinos become degenerate in mass, in which case, one must resum the self-energy 
corrections which lead to the so-called resonant leptogenesis. 

To get a non-zero CP asymmetry, one needs to have complex phases in the mass matrices. 
In the model presented in the last section, we have assumed the right-handed neutrino mass 
matrix is real. Now we relax this assumption, and choose the simplest possibility that the a 
and b parameters have phases a and j3, respectively. Thus the new mass matrix becomes, 

3 iQ 0\ 

M m = ( be* O A*. (18) 
1/ 
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We refit the parameters a, b, a and (3 to reproduce the light-neutrino masses and mixing. 
Clearly, there is no unique set of parameters that are capable of reproducing the phenomenol- 
ogy. One particular set of parameters we use here is 

a = 0.00129 , b = 0.00198 , 

a = -1.808 , (3 = -3.210 . (19) 

Thus the masses of three right-handed neutrinos are then 

M l = 2.27 x 10 10 GeV , M 2 = 3.61 x 10 10 GeV , M 3 = 1.85 x 10 13 GeV . (20) 

We see that Mi and M 2 are fairly close to each other, with 5 = (M 2 — Mi) /Mi = 0.59. The 
Yukawa matrix in the basis in which the right-handed neutrino mass matrix is diagonal and 
real looks like 

/4.8 x 10- 6 exp(0.9z) 0.013 \ 

Yl 3 = -0.022 (21) 

\ 0.0046 exp(0.9i) -0.022 exp(1.6z) 0.472,/ 

where the signs of the phases are chosen to reproduce the right sign for the CP asymmetry. 
Plugging the Yukawa matrix and mass ratios, we find the following CP asymmetries, 

ei = -0.92 x 10~ 5 , 

e 2 = -0.24 x 10~ 5 . (22) 

Here we have also shown the CP asymmetry from the second right-handed neutrino because 
its mass is close to the first one and is potentially important for leptogenesis. The result for 
ei exceeds slightly the bound derived by Davidson and Ibarra [2l| because the masses are 
not so hierarchical. 



B. Effective Out-of-Equilibrium Decays 

In our model, Mi is close to M 2 , and M3 is much heavier. Thus, it is a good approximation 
to neglect the CP asymmetries and lepton number generated from the heaviest right-handed 
neutrinos (those with mass M 3 ). However, since 5 = (M 2 — Mi) /Mi is less than 1, one has 
to consider the full decay and washout effects from the two light right-handed neutrinos. 

The efficiency factor can be calculated by solving the Boltzmann equation for the right- 
handed neutrinos and lepton densities. The result depends on the effective mass mj. In the 
present case, we find, 

fhi = 29.1 meV , fh 2 = 406 meV , (23) 

The effective masses determine the so-called decay parameters = fhi/m* where m* = 
16n 5 / 2 v ^v 2 /(3V5M pl ) = 1.08 x 10~ 3 eV. In our case 

Ki = 27.0 , K 2 = 376.2 . (24) 

Since Ki ^> 1, we are in the so-called strong washout region. In this region, the effective 
factor has little dependence on the thermal history of the right-handed neutrinos. One can 
assume for instance that they are not present in the beginning but are produced purely by 
the inverse scattering process. 
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Since the M\ and M 2 are close to each other, one expects that the existence of N 2 will 
strongly modify the washout of Ni. This situation has been discussed recently in Ref 22l |. 
where analytical formulas have been derived from numerical solutions of the Boltzmann 
equations, 

2 



z B {K x + IC£ 0> )-(K 1 + K 
[l + 21n|g|] 



2 

Sir ; 



^2 = 1 - „ 3 ll - au j—- e-TMi+jJ , (26) 

^(^2 + Af-^).(K 2 + Ki 1 ^) 

where = M\/Tb is the inverse temperature at which the washout effects are minimized 
and k 2 is valid when 5 < 1 |3]. Plugging in the parameters, we find, 

k x = 6.8 x KT 3 , k 2 = 1.3 x 1(T 4 . (27) 

Thus, because K 2 3> ifi, one has K\ 3> /«2- Therefore, the number of out-of-equilibrium 
decays from N 2 is more than an order of magnitude smaller. 

Putting everything together, the baryon asymmetry in our model is 



— 0.96 x 1(T 2 K i e i = 6-° x 10~ 10 > ( 28 ) 



1 



which is close to the observation data. 

We have also numerically solved the Boltzmann equations to check the accuracy of the 
above formula, 



dn Ni 
dz 

drir 



-(A + S) (n Ni - n$) , (29) 



j -e Nl Di (n Nl - n^J - e N2 D 2 (n N2 - ri$ a ) - WN C , (30) 

where z = Mi/T, is in unit of n 7 in a co-moving volume. Di is the decay width measured 
in Hubble expansion rate H . 

Dl ~ KlZ ic^y ° 2 ~ K2Z tc 2 ((i + 5)z) ' (31) 

where /Ci^ are modified Bessel functions. S is the scattering rate of the right-handed neu- 
trinos off the Higgs bosons and gauge bosons. The washout rate W depends on the inverse 
decay, right-handed neutrino scattering, and processes involving the right-handed neutrino 
in the intermediate state (the AL = 2 process). Since we are in the strong washout regime, 
only inverse decay itself can bring the right-handed neutrinos to thermal equilibrium, so we 
may neglect scattering and consider only decays and inverse decays, 

W = ^z 3 [K 1 -JC 1 (z) + K 2 -(l + 5) i IC 1 ((l + 5)z)} . (32) 

The final baryon asymmetry can be calculated via 

' dz i<^i e -J?dz"wv') (33) 

dz' 
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FIG. 1: (color online) The efficiency rate for out-of-equilibrium decays of the lightest (solid line) 
and the next-to-lightest (dashed line) right-handed neutrinos as a function of inverse temperature 
z = Mx/T. 



From Eq. (29), we solve the time evolution of the Ni distribution. 

It was pointed out 2^| that the main contribution to Ki comes from a Gaussian-like peak 



of the integrand around zb- Integrands of both K\ and k 2 are depicted in Fig. 1. We find 
there is such a peak at around zb = 8 > z eq (z eq represents the time that distribution of the 
z-th RH neutrino equals the equilibrium distribution and then follows closely to it in strong 
washout regime). The asymmetries generated before thermalization of Ni are negligible, 
since they experiences longer washout. The dashed line, which shows the integrand for k 2 , 
is much smaller than the solid line for K\. So in our model, we can basically neglect the 
effect of N 2 . The numerical results are 

ki = 5.9 x 1(T 3 , k 2 = 1.4 x KT 4 . (34) 

Compared with Eq. (27), we find that the numerical solution and the analytic approximation 
are reasonably close. 



C. Adding Supersymmetry 

In the presence of supersymmetry, for example in minimal supersymmetric standard 
model (MSSM), the above leptogensis calculation must be modified in several ways. 

First, the lepton number density now must include the contribution from the decay of 
right-handed sneutrinos, 

3C(3)^T 3 ° 

n L = — 2_^{Kiti + Ki€i) , (35) 

i=i 

where, Ki and ?j are the efficient factor and decay asymmetries of the sneutrinos. Because 
of supersymmetry, the second term in the sum is the same as the first term. 

Second, when the lepton number is converted into the baryon number through the 
sphaleron process, one has now n B = —8/23n L . 

Third, the entropy density s = ^g*7r 2 T 3 now has an effective factor 228.75 in MSSM. 
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Combining the above, one has 
344 



Vb 



77165 



Y^iWi + = - 4 - 4 6 x 10~ 3 J2(eiKi + eft) , (36) 



which has a coefficient roughly a factor of 2 smaller. 

Now consider the decay asymmetry. Equation (|16|) for the right-handed neutrino still ap- 
plies, except now we have to take into account the sneutrino intermediate state contribution, 
the decay function becomes, 



F{x) = —\fx 



. l + x 
in 



x — i x 



(37) 



In the limit of large-x, the above becomes 3/ \fx. Thus the asymmetry is a factor of 2 bigger. 
In our model, this is roughly the case: 

e a = -1.78 x 1(T 5 , e 2 = -4.9 x 1(T 6 . (38) 

The decay asymmetry for the sneutrinos <Tj are the same as because of supersymmetry. 

To get the efficiency factors, we turn to the Boltzmann equations. Not only do we now 
have an equation for the right-handed neutrinos, but also for the sneutrinos. Ignoring the 
effects from scattering and taking into account the decay and inverse decay, we have 

drift. n / eq\ 

— - -A (n R - n£) , 



dz ~ lVN > 

-e Ni Di {n Ni - n^J - e^.A (n R - n^J - Wn c , (39) 



dnc 
dz 



where stands for the density of sneutrinos. The above equations are similar to the non- 
super symmetric case, except for the additional contribution from the sneutrinos. Because 
of supersymmetry, the latter contribution to the lepton density is the same as that of the 
neutrino. However, the decay width of the particles is enhanced by a factor of 2, which leads 
to a factor of 2 larger K-factors, 

K x = 53.9, K 2 = 752.5 . (40) 

We numerically solve the Boltzmann equations and find that 

ki = 3.1 x 10~ 3 , k 2 = 4.8 x 10 -5 . (41) 

So Ki are about half of what we find in the non-supersymmetric case. 
Tallying all the changes above, we find the final baryon asymmetry 

T) B = 4.9 x 10~ 10 , (42) 

which is slightly smaller than that of the non-supersymmetric case. But the difference is 
well within theoretical uncertainties. 

In this model, the lightest right-handed neutrino can be produced thermally if the temper- 
ature of the universe after the inflaton decays is higher than 10 10 GeV. However, if one takes 
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FIG. 2: The predictions of Sep, <pi, and 02 plotted against sin 2 6 ) i3 The points are chosen accord- 
ing to the requirement of producing enough leptogenesis and satisfying the 3<r range of nuetrino 
oscillation data as described in the context. 



the cosmological gravitino problem seriously, there is an upper bound on the reheating tem- 
perature Tr < 10 6 — 10 8 GeV when the gravitino mass is in the range 100 GeV < m.3/2 < 1 
TeV [23]. In this case the right-handed neutrinos will be produced non-thermally from the 
decay of the inflaton. For example, one may consider the superpotential 24] 



W^n = ^m^ 2 + g<j)NN , (43) 

where <fi is the inflaton field with mass > 2M\ and g is a dimensionless coupling. The 
reheating temperature is given by 

T R ~ \g\y/wjM*. (44) 

Assuming that the branching ratio of the inflaton decay is of order 1, the produced baryon 
asymmetry is given by 

in 10 e T R 5 x 10 10 GeV 

" * 10 10«GeV ' < 45) 

Taking the reheating temperature Tr ~ 10 7 GeV and ~ 2M\ one can still obtain the 
desired baryon asymmetry of the universe. 

Having shown that our model can generate enough CP violation at high energy for lep- 
togenesis, it is interesting to calculate the size of the CP violation at low energy. The low 
energy CP violation is encoded into one Dirac CP phase 5c p, which is multiplied by sin6 ) 1 3 
in the standard convention, and two Majorana phases <pi and <p2, which appear in the form 
diag(e^ 1 , e^ 1 , 1) in the PMNS matrix. It has been shown by Branco, Morozumi, Nobre 
and Rebelo in |25| that there is no model-independent relation between the CP violation 
at high and low energies. However, we do have predictions of the low energy CP phases 
from our model, and it turns out that they are all small. The scatter plot of these CP 
phases versus the sin 2 6 ) i3 are shown in Fig. 2. Those points are chosen to satisfy the 3a 
range of neutrino oscillation data: 0.7 < sm 2 26>i 2 < 0.92; sin 2 20 2 3 > 0.87; sin 2 ^ < 0.051; 
7.1 x 10~ 5 eV 2 < Am 2 12 < 8.9 x 10~ 5 eV 2 ; 1.4 x 10" 3 eV 2 < Am 2 23 < 3.3 x 10- 3 eV 2 . As shown 
in these scatter plots, 5c p is constrained to be around 3°, and the <fii and $2 are constrained 
to be within 3 degree and 5 degree deviation from —180° and 90°, respectively, indicating 
small CP violations at low energy. The prediction of sin 2 26 ) i3 is shown to be within the 
range 0.06 < sin 2 2^ 13 < 0.085. 
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IV. LEPTOGENESIS IN OTHER SO(10) MODELS 



There are a number of other realistic SO(10) models in the literature which fit well the 
quark and charged-lepton properties, and are consistent with the recent experimental data 
on the neutrino mass differences and mixings. However, they provide very different pictures 
on leptogensis. The key parameter which controls the main features of the leptogenesis is 
the effective mass rh\. For a small fhi, the efficiency factor is large, and one only needs a 
moderate value of the decay asymmetry e to accomplish leptogenesis. For a large values 
of fhi, the out-of-equilibrium decays are rare, and a successful leptogenesis requires a large 
decay asymmetry, which is possible when the masses of right-handed neutrinos become 
degenerate (resonant leptogenesis). 

In this section, we compare leptogenesis scenarios in different SO(10) models and com- 
ment on their strong and weak points. 



AB Model 




The AB model 3| utilizes Higgs fields 10#, 16#, 16# and 45#. The Dirac mass matrices 
of fermions are as follows: 

5 S'e^ 
M u = | -e/3 | , M d = | 5 -e/3 | A D 

8'e i(t > a + e/3 1 

V $n S' n \ f 5 5'e« 

M UD = \ 5 N e Ac/ , Mi = j 5 a + e | A D . (46) 
8' N -e 1 J \5'e^ -e 1 

Notice that the mass matrices Mi and Md are lopsided, producing a large mixing for the 
right-handed down quarks as well as the left-handed charged leptons. Type-one seesaw 
mechanism is used to generate the light left-handed neutrino masses. Since the large-mixing 
angle MSW solar neutrino solution is preferred by experiments, the right-handed Majorana 
mass is constrained to 3 




M VR = -be V e 2 -e \ A R , (47) 



where rj and e are the same parameters as those in the Dirac mass matrices, and a, b and c 
are additional parameters of order 1. 

A set of parameters which reproduce the quark and charged-lepton spectra and mixings 
are, 

e = 0.147 , r] = 6 x 1(T 6 , 



5 = 0.00946 , 5' = 0.00827 , 
a = 1.83 , = 2tt/3 , 
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mu = 113 GeV , mo = 1 GeV 



Given the above, additional parameters, a, b, c and A#, can easily be found to fit the neutrino 
mass differences and mixings. However, the model usually generates a very large m i; which 
in turn produces a very large decay width for the lightest right-handed neutrino. As a 
consequence, the efficiency factor k is too small. To enhance the lepton number production, 
the masses of the two lightest right-handed neutrinos are forced to a near degeneracy, yielding 
a large resonant decay asymmetry. 

In a recent publication, a very extensive search in the parameter space was conducted 
to find a viable leptogenesis in the model ^2$. One of the solutions is described by the 
following parameters, 

7?=l.lxl(r 5 , 5 N = -1.0 x 1(T 5 , S' N = -1.5 x 1(T 5 , 
A R = 2.85 x 10 14 GeV , 

a = c = 0.5828i, b = 1.7670i . (48) 

These parameters lead to the following right-handed neutrino masses, 

Mi ~ M 2 = 5.40 x 10 8 GeV , M 3 = 2.91 x 10 14 GeV . (49) 

The rjB we calculate from these parameters, however, is 2.6 x 10 -6 , roughly a factor of 2 
smaller than that quoted in Ref. |26|. The difference comes from the CP asymmetry of 
the decay. When the masses of the two right-handed neutrinos are close, one cannot use 
the one-loop result in Eq. ()16|) directly. One has to resum the self-energy correction ji?} to 
arrive at 

„ im[(rV%] 



f2 



\m[{Y*Y>y 2l } r N 
8n(Y^Y>) 22 r% + /8tt]* 



where Y' is the vlHNr Yukawa couplings in the mass eigenstate basis of right-handed 
neutrinos, and r N = (Mf - M%)/{M X M 2 ) = -25 is the de generacy parameter. 

It is worth pointing out that although the CP asymmetry tends to be enhanced due to 
the resonance in the case of two lightest right-handed neutrinos being quasi-degenerate, the 
washout effect is also enlarged in this case. Fortunately, in the present model, rh 2 ~ fhi, so 
the effect is not particularly large. The modified numerical results are listed in Table 1. 



B. The BPW Model 

In a model proposed by Babu, Pati, and Wilczek pX| . the fermion Dirac and Majorana 
mass matrices have the following form, 



e + a 1 Au , 1 li 





( \ 


e' 













-e + cr 
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BPW 


GMN 


JLM 


DMM 


AB 


Mi(GeV) 


10 iU 


10 id 


3.77 x 10 1U 


lO 1,3 


5.4 x 10 8 


— e 


2.0 x 10 sin 20 


1.94 x 10 


1.0 x 10 


10 sm 20 


9.4 x 10 


fhi(eV) 


0.003 


0.006 


0.026 


0.1-0.4 


5.4 


K 


6 x 10~ 2 


1.2 x 10" 2 


6.3 x 10~ 3 


10~ 3 


1.4 x 10~ 5 


VB 


12 x lO- lo sin20 


4.97 x 10" 10 


6.2 x 10~ 10 


10~ 9 sin 20 


2.6 x 10~ 10 


sin 2 26»i3 


< 0.1 


0.12 


0.06-0.085 


0.014 - 0.048 


0.01 



TABLE I: Predicted mass M\ of the lightest right-handed neutrino, CP asymmetry e, effective mass 
fhl> efficiency factor k and baryon asymmetry t]b, and #13 in various SO(10) models. The order 
is arranged according to the size of rh\. 



-3e' \ / -3e' + rf 

M„ D = I 3e' -3e + a \a v , Mi = 3e' + T]' -3e + r] \ \ L > . 
3e + a 1 / \ 3e + 77 1 

fx z\ 

= y \A R . (51) 
\z y 1/ 

One set of parameters which produces good phenomenology without CP violation is 



(7 = 0.110, 77 = 0.151, e= -0.095 



/ = 4.4 x 10~ 3 , e' = 2 x 10" 4 , 



a; = 10- 4 , y=-l/17, z = 1/200 



A l/ (M X ) = 120 GeV , A D (M X ) = 1.5 GeV , A R = 10 15 GeV . 

With the above, a number of successful predictions follow, including the masses for bottom 
and down quarks, CKM martix elements, and the atmospheric neutrino oscillation parame- 
ters. 

However, the solar-neutrino mixing angle from this model comes out too small: sin #12 = 
0.04, in contradiction with the experimentally-preferred large-mixing angle solution. To 
remedy this, a small intrinsic mass for lefthanded neutrinos is introduced, which stems from 
operators k^IGi^IG^IGf^OhIOh/Mq. With this modification, the solar mixing angle 
changes to sin 2 2#i 2 — 0.6 |28| |. 

The result of leptogenesis in this model is summarized under "BPW" in Table I. Instead 
of showing a range of results, as in Ref. j2^|, we take mi = 3 meV, and the efficiency factor 
is then about 6 x 10 -2 . The baryon asymmetry is about 12 sin 20 x 10~ 10 , where is a 
CP- violation phase. With a reasonable choice of 0, the experimental tjb is produced. 

Clearly the model does not provide a tight constraint on the relation between the low- 
energy neutrino properties and leptogenesis. Introducing an intrinsic mass term for the 
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light neutrinos might be physically motivated. However, it dilutes the relation between Mr 
and low-energy neutrino observables. In an extreme case, the constraint on the right-handed 
neutrino properties will be lost if the low-energy neutrino mass matrix is entirely "intrinsic" , 
i.e., of non-see-saw origin. A mild reflection of this "decoupling" is that rh\ in this model 
is particularly small, which is hard to achieve in a complete first-type see-saw model which 
fits the low-energy data. A small rh\ yields a large efficiency factor which certainly aids the 
leptogenesis here. 



C. The Minimal 126-Higgs Model 

This model 0, H3| (referred to as GMN model in the table) uses the 126// to break B — L 
symmetry and the 10//, 210// to break gauge symmetries [3], and allows all couplings to be 
complex. 10// and 126// are used to give fermion masses through superpotential 

W Y = hjiij^jHw + fij^jAjsQ , (52) 

So the following mass matrices are obtained: 

M u = h + f, M d = hn + fr 2 , 

Mi = Tin - 3/r 2 , M UD = h- 3f , (53) 
where h and / are matrices 

h = h* cos a u sin (3 , / = f*e iryu sin a u sin (3 , 

and the parameters are 

r 1 = ^cot/3, r 2 = e^--)^cot/3. (54) 
cos a u sin a u 

Hence, we can get a sumrule for quark and charged-lepton masses 

k M, =r M i + M^ 
m T rrib m t 

The seesaw formula in this model is given by M Vh = fv\ — M VD {M VR )~ X {M VD ) T and 
the calculations are done assuming that the first term dominates over the second. The 
left-handed neutrino Majorana mass matrix satisfies a sumrule 

M VL = a(M t - M d )* , (56) 

The right-handed Majorana neutrino mass matrix is given by choosing the B — L (seesaw) 
scale vb-l = 2 x 10 14 GeV and ^ u = 0, sina u ~ sin/5 ~ 0(1), 



e 



M VR = fvB-L = f* . - ■ ^ b-l ^ (2 x 10 14 GeV)f • (57) 
sin a u sin p 

Therefore, the left- and right-handed neutrino Majorana masses have the same texture, and 
there are some nontrivial relations between the low energy phenomena and leptogenesis, 
when the first term is assumed to dominate in M VL . 
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We use the following parameters and matrices, 

k = -0.846 , r = -1.846 , 



(58) 



ri = 0.0116 , r 2 = 0.00337 , (59) 



/ 0.1768 -0.0207 + 0.0038i -0.058 + 0.097i 

/••-•>',, — .1/,; / , 

h=^-^ -= -0.0207 - 0.0038i 3.49 -1.332 , (60) 



r 2 - ri 



-0.058 + 0.097i -1.332 94.7 



0.165 0.0714 - 0.0133z 0.2 - 0.335i 

f = M u ~h= | 0.0714 + 0.0133i -3.161 4.596 J . (61) 

0.2 + 0.335i 4.596 -12.28 

The predicted baryon asymmetry through thermal leptogenesis is listed in Table 1 as 
"GMN" . 

This model is also characterized by a small mi, which is not constrained by the low-energy 
neutrino mass spectra because the first type of see-saw mass contribution is assumed to be 
small. Note that the neutrino Dirac mass texture in this model is completely untested at 
low energy and can only be effective in the right-handed neutrino decay process. 

D. The DMM Model 

The DMM model |9( is an extension of the GMN model by enlarging the Higgs sector to 
include 120(A) which gives an additional contribution to the fermion mass matrices through 
the coupling 

^120 = iK^D , (62) 

where h'^ is an antisymmetric matrix due to the S'O(IO) symmetry. Furthermore, the 
10 and 126 couplings are chosen real and the 120 imaginary by using a Z 2 symmetry. 
In this case there are six pairs of Higgs doublets: ifd = (H^°, A d , A\, Ad, Ad, $d), (p u = 
(H™, D\, D\, A u , A u , $„), where superscripts 1, 2 of A u ^ d stand for SU(4) singlet and adjoint 
pieces under the G422 = SU(4)xSU(2)xSU(2) decomposition. 
The MSSM Higgs doublets are given by 

H d = UlMa , (63) 
H u = V* a ((j) u )a , 

where a — 1, ..6, U and V are unitary matrices which diagonalize the Higgs mass matrix. 
The Yukawa coupling matrices for fermions are given by 

Y u = h + rj+riti , (64) 

Yd = nCh + f + h') , (65) 

Y e = ri (h-3f + cji') , (66) 

Y v = h-3r 2 f + c u h' , (67) 
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where the subscripts u, d, e, v denote for up-type quark, down-type quark, charged-lepton, 
and Dirac neutrino Yukawa couplings, respectively, and 

h = V n h, f = U u /(V3n)f , h' = (u 12 + u 13 /Vs)/nh' , (68) 
U n V 15 V 12 - V 13 /V3 

n = Wr r2 = n W.> r3 = ri u 12 + u 13 /V3> m 

U 12 - V3U 13 V 12 + V3V 13 



n 1Z v x t- . (70) 



6 U12 + U 13 /V3' " U 12 + U 13 /V3' 

The light neutrino mass is obtained as 

m^ ht = M L - M°M^{M°) T , (71) 

where = Y V (H U ), M L = 2V2f(A L ), and M R = 2y/2f(K R ). 
It has been shown [9| that the Dirac neutrino coupling matrix 

/ 0.002 0.003 exp(-1.54i) 0.0026 exp(-0.344i)\ 
% = j -0.0167 0.021 exp(-1.53i) 0.025 exp(3.37i) , (72) 
\ -0.229 0.417 exp(4.70i) 0.422 exp(-3.019 i) J 

gives a good fit to Am^ tm , Am s 2 | and mixing angles with the lightest right-handed neutrino 
mass Mi ~ 10 13 GeV. In this model there is a correlation between U e3 and V u b as well as U e3 
and Aml ol /Am\. The former imposes upper bound on U e3 , while the latter gives a lower 
bound. The predicted values of U e3 , lepton asymmetry and the washout factor are given in 
Table 1. 



V. CONCLUSION 

In summary, we have discussed leptogenesis in a SUSY SO(10) GUT model for the fermion 
masses and mixings, which is developed from the original lopsided model of Albright and 
Barr We have done a detailed analysis of the washout factor and find that the model 
predicts a value for the baryon to photon ratio (tjb) of the universe in good agreement with 
observations from WMAP as well as the requirement of a successful nucleosynthesis. We 
then compare with the same predictions for other successful SO(10) models in the literature. 
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DE-FG02-93ER-40762 and by National Natural Science Foundation of China (NSFC). R. 
Mohapatra and S.Nasri are supported by National Science Foundation (NSF) Grant No. 
PHY-0354401. 
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